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Abstract We study different types of radiative decays involving /o(980) and ao(980) mesons within a unified 
ChPT-based approach at one-loop level. Light scalar resonances which are seen in tttt, tti], KK channels of 
(^(1020) radiative decays and in J/tp decays are responsible for key questions of low-energy dynamics in the 
strong interaction sector, and decays <^(1020) ^ 7ao(980), <?!>(1020) 7/o(980), ao(980) ^ 77, /o(980) 
77 are of interest for current experimental programs in Jiilich, Prascati and Novosibirsk. From theoretical 
point of view it is important to verify whether light scalar mesons are members of some flavor octet or 
nonet. We find a value of mixing angle dictated by consistency with experiment and coupling structures of 
ChPT Lagrangian. Decay widths /o(980)/ao(980) 7p(770)/lj(782), which are not studied experimentally 
yet, are predicted. We also obtain several relations between widths, which hold independently of coupling 
constants and represent a fingerprint of the model. 



X 



PACS. 11.30.Hv Flavor symmetries 
Properties of mesons 



12.39.Fe Chiral Lagrangians - 13. 30. Eg Hadronic decays - 14.40.-n 



1 Introduction 

The scalar mesons ao(980) (/"^(J^^) = 1^(0++)) and 
/o(980) {I^(JP^) = 0+(0++)) have been discussed for 
more than 30 years. The shape of tttt (and irrj) invari- 
ant mass distribution in different reactions points to these 
resonances. The promising source of information on scalar 
mesons are radiative decays in which scalar mesons are in- 
volved. Much experimental attention has been paid so far 
to processes 0(1020) 700 [I] and 0(1020) 7/0 [2J 
(see also [3l|4|) due to motivation put forward in Re- 
cent example of a model describing such features in the 
rare </> — 7^ jirr] {jittt) decays is chiral approach 
with derivative coupHngs [6j. Among other well-known 
processes involving scalar resonances one can think of 
J/ip 0/0(980) ^ 07r7r (and ^ (l)KK) studied at BES [7] 
and nucleon-nucleon (as well as deuteron-deuteron) reac- 
tions with various hadronic final states. The transitions 
ao 77 and /o 77 are relevant for numerous re- 
actions, where two-photon interactions produce miscel- 
laneous hadronic final states. Many experiments involv- 
ing 77 TTTT (or mj) have been carried out or are being 
planned. 

The properties of scalar mesons are not well un- 
derstood. Nevertheless the dominant decay channels are 
known to be tttt for fo meson and ttt] for ao, and the to- 
tal widths are in between. 40 and 100 MeV. The decays to 
strange mesons ao KK and fa KK are dynamically 
allowed, though the masses of ao(980) and /o(980) may lie 



sHghtly below the KK threshold. The masses of ao( 
and /o(980) are approximately equal. 

The internal structure of light scalar mesons is also not 
clear. Recent advances in understanding of their struc- 
ture are presented in review [8]. Most of studies show 
that light scalar meson structure can not be explained 
in simple quark models. This is probably related to a spe- 
cial role played by these mesons in low-energy dynam- 
ics of strong interaction [OlfTO]. Namely, scalar fields can 
be viewed as the Higgs sector of strong interaction, i.e. 
their non-zero vacuum expectation value leads to chiral 
symmetry breaking and directly refiect the structure of 
quark condensate in Quantum Chromodynarnics (QCD). 
Some authors emphasize proximity of the KK threshold 
to ao and fo masses that favors presence of the molecular 
KK componejit [Tl] (for recent calculations implementing 
molecular KK model see [T2 lll3pi4j .) 

Another approach to light scalar meson features 
is unitarized ChPT [T^. There the inverse amplitude 
method_|16| was employed to describe elastic tttt, tttj, Ktj 
and KK scattering data. The radiative decays in question 
were evaluated through final state interaction of scattered 
particles. The unusual large- iVc behavior of scalar reso- 
nances was recently summarized in [17] (see also references 
therein and original paper [T8]l. 

The decays /o/ao(980) 7 p(770)/w(782) are similar 
to decays ^ 7 ao (7 fo)- The interest to these processes 
was initiated in [19]. Apparently they can be explained 
in terms of the same matrix element (describing process 
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with vector (V) and scalar (S) particles, and photon in 
initial/final state) with SU{3) flavor modiflcations reflect- 
ing the type of vector particles. Decays S may be 
studied experimentally in Jiilich (with ANKE and WASA 
at COSY) |20] and possibly in Frascati [21] (with KLOE 
at DA^NE or its upgrade) and at BES. 

Various phenomenological models [F'.'6'.'T2 |I131I14P22U231 
[24] have been applied to calculation of these decays. At the 
same time a consistent description in framework of Chi- 
ral Perturbation Theory (ChPT) with pre-existing vector 
and axial-vector mesons p5] is lacking. This is an effec- 
tive theory of the strong and electromagnetic interactions 
at energies below 1 GeV and has symmetries of the un- 
derlying QCD. Strictly speaking ChPT is an expansion in 
series p'^/A^, m^/A^, where p is momentum, m is mass of 
pseudoscalar mesons, and chiral symmetry breaking scale 

is of order 1 GeV. Thus, formally, the range of energies 
for scalar mesons is on the border of ChPT applicability. 
Nevertheless, it is clear that suitable effective Lagrangian 
for scalar mesons has to have much in common with ChPT 
Lagrangian, because the coupling structures are guided by 
the chiral symmetry. There are many successful applica- 
tions of this theory at energies about 1 GeV that make a 
useful background for employing it in present problem. 

In general, ChPT does not specify internal structure 
of interacting particles. The model [25] only assumes that 
the scalar flelds belong to SU{3) flavor octet and singlet. 
This Lagrangian is written down in Appendices [A] and [Bl 
in particular, describes interaction of pseudoscalar and 
vector mesons, and - interaction of scalar mesons with 
pseudoscalars. We test the singlet-octet mixing scheme for 
the lightest scalar meson nonet 

{ao ^Ss, 
/o =5"™s C0S61-58 sin6', (1) 
a =5«'"s sin 6* 5*8 cos 6*, 

where 5*3 is the neutral isospin-one, Ss is isospin-zero 
members of flavor octet and S"'™^ is flavor singlet. 6 is 
the octet-singlet mixing angle, and a — /o(600). In par- 
ticular we are interested in whether ao(980) and /o(980) 
suit for members of this nonet. In principle this may not 
be the case (see for example argumentation in [26]) and 
therefore should be verifled. Radiative decays may help to 
clarify this important issue. 

The present paper considers decays S 77, 
0(1020) — > -fS and S ^V. We suppose that the under- 
lying dynamics of all above decays has much in common, 
namely that the loops with pseudoscalar mesons form the 
dominant mechanism. This assumption is consistently im- 
plemented in Lagrangian |25j and Section[2]presents calcu- 
lation of decay amplitudes. On this way we prove cancela- 
tion of divergences and gauge invariance of the amplitudes. 
Along the calculations we use the dimensional regulariza- 
tion method, see Appendix [D] for a brief overview of the 
method and list of basic formulae. Some details of loop 
integrals calculation and their analysis are also presented 
in Section [21 

The are six coupling constants in Lagrangian {Fy , Gy , 
Cd, Cm, Cd, Cm), and estimation of their values is carried 



out in Section [31 Under assumption of the resonance sat- 
uration the coupling constants may be expressed in terms 
of chiral LEC's [25]. Available experimental data provides 
certain constraints on these couplings. 

After fixing the parameters we calculate the widths 
of various decays with light scalar mesons and compare 
them with available data and predictions of other models 
(Section [3|). We compare pion and kaon loop contribu- 
tions to the decays with /o meson in initial/final state 
and demonstrate the importance of pion loops in fa jp 
and /o 77 decays. 

The virtual photon case, which is important for fur- 
ther applications of the present model, is outlined in Ap- 
pendix [Cl 

2 Formalism for radiative decays amplitudes 

2.1 One-loop diagrams and chiral counting 

From Lagrangian terms IjA.ip and l|B.4p one obtains the 
sets of one-loop diagrams shown in Figs. [H [2l and [3l In 
the present approach we have no tree-level diagrams for 
the radiative processes. Therefore the lowest-order ampli- 
tudes consist of one- loop diagrams. The corresponding set 
of diagrams for ao//o 77 decay with pseudoscalar me- 
son in the loop is shown in Fig. [H This set of diagrams is 
complete since it is obtained from Lagrangian which car- 
ries the chiral power not less than the chiral power of any 
diagram. 

The following rules [27] are used to count the chiral 
power of any diagram. These counting rules provide one 
with a guiding idea of which diagrams should be included 
and which should not, when forming a set of relevant di- 
agrams at any given order. Pseudoscalar fields scalar 
fields S and vector fields (in tensor representation) p^i/, 
w^i^, ^/if carry zero chiral power 0{p'^); derivative or ex- 
ternal source (like electromagnetic field i?^) has unit chi- 
ral power 0{p); pseudoscalar-meson mass (m^, uik) also 
carries unit power 0{p) (so that the mass matrix x in Ap- 
pendix[A]isO(p2)). The propagator of pseudoscalar meson 
is counted as 0{p~'^). 

All coupling constants in Lagrangian (cd, Cm, Cd, Cd, 
Fv, Gy) are 0(j3°), and the power of any vertex is deter- 
mined only by the structure of the corresponding term in 
Lagrangian. In addition, the loop integration over parti- 
cle momentum adds 0(p*). Applying these rules one can 
show that each diagram in Fig. [H [2 and [31 has chiral power 
0{p'). 

2.2 Radiative decays ao/ fo "f"f 

Consider the amplitude of ao/ fo 77 decay. Let scalar 
meson have 4-momentum p, photons have polarization 
vectors ej^^'' and e;y \ and 4-momenta and Suppose 
that a positive charge runs clockwise in the loop. First, we 
write down amplitude for the first three diagrams (a-c), 
as this set of diagrams is often used in various approaches 
for radiative decays [SllGllTO]. 
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Figure 1. Diagrams for decay of scalar meson (dotted line) into two photons (wavy lines). Pseudoscalar mesons (dashed lines) 
in the loops are: {K^K~) for ao(980) decays, and (K^K^), (tv^tt^) for /o(980) decays. Solid (non-derivative coupling) and 
crossed (derivative coupling) blobs represent 0{p^) vertices. Arrows mark the places, where form factors of pion and kaon would 
arise for virtual photons (see discussion in Appendix [C|l . 



The invariant amplitude Mate corresponding to di- 
agrams (a-c) in Fig. [1] is expressed through the tensor 

rpfJ.L' 

S(P)77- 



„M . - e(i)* e^^'* T'"' 



(2) 



The indices S{P) indicate that the scalar meson of type S 
decays in two photons via loop consisting of two interme- 
diate pseudoscalar mesons of type P with mass mp, i.e. 
S{P) = {ao{KK),fo{KK)Jo{n7T)}. 
Explicitly we find 



e gspp 



+r^(gi,g2)+/^"^fe,gi)}, (3) 
'^'^ {2l-p-qiY{2l~qiY 



(2^)4 

X lAlAl^pAl^q 



where Ai 



. For gspp^s we refer to La- 



grangian (|B.4p in Appendix [BJ Changing the integration 
variable V = p — l (and V = p + qi — I) and assuming that 
the possible divergence of the integrals is not higher than 
the logarithmic, one can prove the gauge invariance of the 
amplitude: 



rjiUU 



S(P)77 



rr-ifll' r\ 



(4) 



Making use of the change I' = p — I we deduce a useful 
relation 



r^{qi,q2) = I^''{q2,qi) 



(5) 



which is connected with Bose-symmetry of the final pho- 
tons. 

By means of the Feynman parametrization and dimen- 
sional regularization method (see Appendix[D]) the expres- 
sion |(3]) is reduced to 

T^fr„..., = -H!£^|g''- f\xHm%-p\il - x)] 



-SiPh^ - (4^)2 f2 ^9 



'^9^" J J xdxdy\n[C{x,y;qi,q2)] 

\dxdv^^^^^^^^^^^^^\ 
C{x,y;qi,q2) > 



(6) 



A'^{x,y-quq2) = (gi[2x(l - y) - 2] + 2xy p - 92)"' 
X {q^\2x(\-y)-\\^2xypY, 
C(x,y; 91,52) = x(x - 1)(1 - y) + p^ xy{x - 1) 

-ql -y) + ^n^- 

The divergent parts of diagrams (a-c) in Fig. [1] cancel each 
other. 

For real photons in question 

ql = ql = 0, 
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Table 1. Values of the loop integrals. Assumed physical values 
of scalar meson masses: M/^ = 980 MeV and Mao = 984.7 MeV 





= 0.0749 + 0.2441 








= 0.5510 - 0.2441 




= -0.63 




= 0.5507- 0.246 4 


a^(mi;M? ;0) 
^ {ml; Mho) 
^{ml;Mh,0) 
^{ml;Ml-,Ml) 
^{ml;MX;Ml) 


= -0.57 


= 0.3545 - 0.56641 


= -0.0529 


= 0.128 - 0.01691 


= -0.4048 


^{mj,;Ml;Mi) 


= -0.3988 


nml;Ml;Ml) 


= -0.3466 


nml;Ml;Ml) 


= -0.3407 



These conditions simplify equation ^ to 

-4ie^ gspp 



Here we define 



[9 



xdxdy 



QUI 
91-92 



1 



)'F(ml;p^;0). (7) 



(8) 



p'^xy{x — 1) 
1 mp dx , , , , »■ 



Jo x-l 

1 Tod dx 



2 



, o T , 1 Tod f dx 



In 



l + .x(a;-l)£^ 



TO, 



1 



1_4I 



1- a/1 -41 



/to tf'(m2;/;0) = TT^ln 



/to i/'(TO|-;p2;0) = 0. 

The integrals ^{nip;p'^; 0) are calculated numerically and 
presented in Table [H The scalar-meson invariant mass 
^^2-^1/2 |g equal to Mg - the mass of /o, oq (and a for com- 
pleteness), while TOP is equal to the mass of pseudoscalar 
meson tt or /f in the loopEl- 

Now we consider diagrams (d-i) in Fig. [TJ At first 
glance, these diagrams, which include derivative coupling 
for scalar mesons q are more complicated due to momen- 
tum dependence of the SPP vertex. In fact these diagrams 
can be treated similarly to the previous case. To demon- 
strate this let us define 



I' S(P)77' 



(9) 



^ When working with integrals (H} it is convenient to use the 
identity /g^(l — 2x)f{y)dx — for any function f{y), where 
y = x{l- x). 

^ We also include here diagram (i), though it has no deriva- 
tive coupling in SPP^y vertex. This is convenient due to its 
cancelation with contribution of diagram (f). 



where symbol "hat" is used hereafter to indicate the deriva- 
tive coupling. Next use the identity 

and change integration variables as above in order to com- 
bine six terms in Md-i in such a way that the contribution 
of diagram (i) cancels the contribution of (f), and dia- 
grams (g) , (h) are cancelled by part of (d) and part of (e) . 
In this way the derivative coupling amplitude T^^p^^^ is 
related to the non-derivative coupling amplitude T^^p^^^ 



.(1).(2) fM"^ 



9SPP 
gspp 



(to|,-pV2) e«^(^^r^Ti 



5(P)77- 



Combining contributions of all diagrams in Fig. [T] one 
obtains the total ©(p^) invariant amplitude 



-4ie"^ 



(47r)V^ 

.(!)■ 



^iml;p^;0)AKip^) 



(10) 



.(2)* 



9i 



/o-»77 



qi ■ 92 

[B^{p^)^{ml;p';Q) 

,2^,7>/™2 . „2 

(1) 



+BK{p')1'{mi,;p'M 



(11) 



. ^(2)* 



92 e'^^* • qi 



qi ■ q2 



where 



Ak{p^) = gaKK{m]i -p^/2) +gaKK, 

Bk{jp^) = g/KKiml; - p^/2) + g/KK, 
B^{p'^) = gf-Kniml - p^/2) + gf^^. 

2.3 Radiative decays 0(1020) — * 7 ao//o 



(12) 



Let vector meson 0(1020) with momentum Q decay into 
scalar meson ao(980) (or /o(980)) with momentum p and 
photon with momentum q, i.e. ^(Q) 7(9) + ao/ foip)- 
Diagrams corresponding to these reactions are shown in 
Fig. [2 

Let polarization vector for the (/)-meson be E^, and 
that for the photon e^. Apparently q'^e^ = 0, Q'^E^ = 
0. We describe vector meson 0(1020) by antisymmetric 
tensor field carrying the indices fiX. Thus we employ the 
normalization for one-particle matrix element [25] 



(0 I 0^a(O) I 0, Q) ^tM^'iQ^Ex - QxE^] 
The invariant amplitude reads: 

QfiEx — QxEfj, 



(13) 



iM^^jS = C ^- 



75 ^<P^-fS 



(14) 
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(A) ^ 



(B) 



(E) 




(C) 



(F) 



(D) 



(G) 



Figure 2. Diagrams for </>(1020) to 7ao(980) or 7/o(980) decays. Pseudoscalar mesons here are: [K'^K ). Solid and crossed 
blobs stand for 0{p^) vertices, the latter indicate derivative coupling terms. 



where tensor with "hat" for diagrams (A-D) is 



= T 



T 



T 



leGvgsKK 



T, 



c 



T; 



J (2^)4 
X {2l-QyAiAi 
-legSKK f dH 



{Q^l^ - - A)) 



{Gv + ^{Fv ^2Gv)q^)g'' 



(27r)4 
1 



in^X) l{l-Q + q)AiAi-Q+. 



9! 



^Xv 



D 

e-GvgsKK 



dH 



{Q^^ - (/i ^ A)) 



J (2^)4 

X {21 ~ qY{l -q){l~ Q) AiAi-^Ai-Q, (15) 
and tensor without "hat" for diagrams (E-G) reads 



rp^Xu rp^Xu 



fiXiy 



dH 



V2/4 J (27r) 



'^Xu 



'legsKK 



g'^Gv Q' 



+ -{Fv~2Gv)q^)-{f^-^X) 



fiXiy 



fiXv 



-eGvgsKK f dH 



V2/4 J (27r)4 



AiAi-gAi-Q{2l- qY 
(16) 



The consideration shows that divergent parts of the 
ampUtudes lfT5|) and (fT6| which do not cancel are pro- 
portional to (Fv - 2 Gv)ig''^q^ ~ g^^q^)- Therefore we 
employ the relation 



Fv — 2 G\ 



(17) 



between electromagnetic and strong coupling constants of 
vector mesons (see Appendix [A]) in order to make the 
amplitudes finite. Actually this relation does not follow 



from the chiral symmetry. However it naturally appears 
in alternative approaches, Hidden Local Gauge Symme- 
try Model IMi and massive Yang-Mills theory [29\. This 
aspect has been addressed in [30]. In Appendix [A] we dis- 
cuss accuracy of (flTl) based on experiment. 

Making use of (fl7| and identities 



i{Ai^q + Ai^Q) = Ai^qAi^Q\p^ -2rnj, + 2{l ^ q){l ^ Q)] 
i{Ai + Ai^Q+g) = AiAi^Q+g [p^ - 2ml + '^1(1 - Q + q)] 



in (|T5|) one can prove that 



rpuXiJ _ gsKK 



gSKK 



xe., i- 



P \rpflXv 

QfiEx — QxEfj 



= 0. (18) 



Finally 



iM 



gSKK , 2 
gSKK 



EL 

2 



(19) 



In calculation of the amplitude the Feynman 
parametrization and dimensional regularization method 
are applied (Appendix [D| . Then eq. ifTOj) . with the use of 
(fTBl). reads 



-iM 



-leGygsKKQ^ 
a/2/4(4^)2M^ 



(20) 



gSKK ( 2 

(m^ 



)) 



gSKK 

X \{Q-e*){q-E)h^{e* ■E){h-2h) 
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where 

12 = c?x ln(m|. — p^a;(l — a;)), 

Jo 

13 = J J xdxdy\n{ml^ — Q^xil — x) -\- 2xy{l — x)Q ■ q), 
»i 



/i ^ / / xdxdy _ 



xy{l - x) 







mj^ — Q'^x{l — x) + 2xy{l — x)Q ■ q 



AQ-q 



{h-2h), 



h - 2/, = 1 - 



X In ■ 



^^ '^K-Mlxil-x) 
''{Ml-p^)x{l~x) 



(21) 



(22) 



•■K 



' p'^x{\ — x) 



mj^ - M^x{l - x) 



2^{ml,p\Ml). 



In terms of tf'(m|-,p^, M|) the invariant amplitude l(20l 
reads 



(23) 



1 



Li ■ q 

^Kip^) for (t)~^KK^ao^ 
BKip") for (j)^KK^foj 



Compare definition ((22l) with that of if'(TOp,p^,0) 
in ([H]). Real and imaginary parts of ^{m'j^,p'^ , M^) at 
p2 — ^2 g^j.g 



Re^iml,MlMl)^ - -- 



1 ml^Mlx{l-x) 
dx- 



2' 2]^ '^''jMf-M^)x{l-x) 



X In 



m|- — M^x{l — x) 



Im^{ml,MlMl) 



mj^ - M^x{l - x) 



'1 



Ml - M2 



Ml 



Ml 
Ml 



1 -4 



hp 



I illK 
Ml 



(24) 



Numerical calculation of W(rnj^, M^ , M^) leads to values 
shown in Table [T] (see also [5] for analytic expression of 
the integral 



the matrix element for these decays is very similar to that 
in One can replace 0(1020) by p(770) (or w(782)), 

take into account fiavor SU{3) factor in the VPP ver- 
tices, and select the pseudoscalar particles in the loops 
allowed by symmetries of the strong interaction. For rele- 
vant SU(2>) relations see AppendixO Taking into account 
that ao(980) and ^(782) do not couple to two pions, one is 
left with K'^K~ loop for ao ^ and /o ^UJ decays, 
and both tt+tt^ and K^K^ loops for /o 7P decay. 
The matrix elements read 



-leMvGv 



X AKip") nmhP^;M^ 



e*-E*^-^{Q-e*){q-E*) 
Q-q 



(25) 



3M„G 



E* 



Q-q 



X {BK{p^)M^{ml-p'-Ml) 
+2B^{p^)^{ml-p^-Ml)), 



{Q-e*){q-E* 



(26) 



iM 



'leM^Gx 



X BKip") ^{ml-p^-Ml), 



E*--^{Q-e*){q-E*) 



(27) 



where notation for momenta and polarization vectors is 
the same as in Section! 



The loop integrals W{ml;p'^; M^) and if'(m|.;p2; M^) 
can also be defined by (|22l) . Their numerical values are 
shown in Table [TJ 

To make a correspondence with results of [22] we can 
write the loop integrals for /o/oo — *■ 7 p/(^ diagrams in 
Fig. [3] in the form 



Re^{ml;M^;M^ 



-b 



In 



2 1 + y/T^I/b 



-In 



2 i + yi^lA 



a ( I 1 + - 4/6 

2{a~h)\^ ' 1-Jl- 4/6 



1 - a/1 - 4/a^ 



(28) 



2.4 Radiative decays /o/ao ^ p/uo 

Decay of scalar meson into vector meson with radiation 
of photon {fo/ao 7 p/uj) in the lowest order is repre- 
sented by O(p^) diagrams shown in Fig. [3l The vertices 
follow from L"^ in (|a7iP and in l|B.4p . The structure of 



1 -I- ^1 - 4/5 



/,.^«;M-M,^)^-^(ln^ 



-71^4/6 



In 



i + v/r^4A 



1-^/1^4/^ 
2(^(v^-^^)' 



(29) 
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(B) 




,>, 

(C) 



..^..4 'p=>= ■■>■■• t 

(E) (F) ^^~^'\ 



(D) 



(G) 



Figure 3. Diagrams for /o(980)/ao(980) — * 7p(770)/lj(782) decays. Pseudoscalar mesons in the loops are: {K'^K~) for all 
decays, and {tt'^iy~ , K") for /o(980) — » 7p(770). Solid and crossed blobs stand for 0{p^) vertices, the latter are derivative 
coupling terms. 



with a = Mp/m^ and b = Mg/m^ for pions in the loop. 
For kaons in the loop one finds 



1 



a — b\ 



, . 2 ■ 2 Va 

arcsm — arcsm — — 



, . v/47&~Tarcsin2^ 

/ — ; 2 ■\/^\ 

— V4/a — 1 arcsm ~^J: 



(30) 



where a = My/m\ and 6 — M^/m^. Numerical values 
obtained from these analytic expressions agree with those 
in Table [U obtained from direct numerical integration. For 
the integrals ([8]) and (|22l) one can also deduce analytic 
expressions from ^22j by choosing appropriate a and 6. 



The invariant mass of scalar meson is 



For decays into KK in ([341) and ([35| one includes 
factor 2, (as KK = K+K~,K°K°), and (1 + 1/2 (2 x 
1/2)2) ^ 3/2 for TTTT in 1 - from charged pions, 1/2 

- from the identity of neutral pions, (1/2)^ - because the 
neutral pions interact two times weaker than the charged 



ones (tt 



27r+7r ), and 2 is the symmetry factor 



in the vertex with two identical neutral pions. 

The widths of our premium interest are built up of 
AkIp'^), Bxip^) and B^{p'^) loop integrals W and 

phase-space factors. Thus, through relations l|12p . they 
depend on Lagrangian couplings Cd, Cm, Cd, Cm and singlet- 
octet mixing angle for scalar mesons (see Appendix [B|) . 

The widths for ao ^ 77, /o — > 77 decays read 



3 Results and discussion 

3.1 Widths and estimates for chiral couplings 

First of all there are direct decays, which can be described 
from l|B.4p at tree-level. They represent the dominant 
channels: oq irri for isotriplet and /o — > tttt for isos- 
inglet scalar mesons 



1 



m2)2 



4p2 



J TT 



1. 1 



rfo^TTTT = (1 + 2) 

here A^,,(p2) is introduced by analogy with (fT2l) 

A-T.-nip'^) = 9a^7j{m'^ + ^1- P^)/2 + 9aTi7j, 
and the decays ao KK and /o KK: 



^V'pV4-<^|Ak(p^)P, 
1 nrr. r i 



(32) 
(33) 

(34) 
(35) 



r, 



1 



00^77 



327rv^8^V^ 

„4 



\AKip^) ^(m2,;p2;0)|2,(36) 



1 



327rV^87r4/.* 

+i?,(p2) ^^(^2.^2.0)12^ 



(37) 



In deriving l(36|) and l(37|) the formula for the width 



Ts^^^ = 1/(2 X 167rAf,)|Xs^77l' 



(38) 



is used, with the amplitude defined in ifTO]) and ifTTj) . and 
symmetry factor 1/2 for two identical photons in the final 
state. Further, the sum over polarization states A of the 
photon is performed using 



A=±l 



(39) 



under condition that polarization vector is contracted with 
the conserved current. 
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Table 2. Particle properties (data from PDG [32]) which are needed in calculations 



Particle 



Mass 
(MeV) 



Width 
(MeV) 



Major hadronic 
decay channels 



ao(980) 
/o(980) 
a = /o(600) 
a = /o(600) 

^± 

0(1020) 



l-(0++) 
0+(0++) 
0+(0++) [32] 
0+(0++) (333 

i-(o-) 

1/2(0-) 

o-(i--) 



984.7 ± 1.2 

980 ± 10 
400 - 1000 
513 ± 32 
139.57018 ± 0.00035 
493.677 ±0.016 
1019 ±0.019 



50 - 100 
40 - 100 
600 - 1000 
335 ± 67 
mean life 2.6 x 10" 
mean life 1.24 x 10" 
4.26 



TTTT 
TTTT 
TTTT 



K+K-, KlK%, 
pn ± Stt 



The widths for (/)(1020) meson decays are 
1 2 



AttM^ 3 



(1 



Ml 



1 



<P~*lfo 



4(1 



AttMj, 3 



P ' 
Ml' 



V2eM^Gi 



(40) 



(41) 



Equations l(40|) and l(4T1) are derived from general expres- 
sion 



r, 



0-*7ao//o 



■fl - — ) 



(42) 



with assumption (fT7| . The amphtude M^i^-ys in given in 
l(23l) . The factor of 2/3 in ^ and {Hj) comes from the 
sum over photon polarizations and average over vector- 
meson polarizations A by means of 



A=0,±1 



(A)^(A)* ^ 



~Mf- 



(43) 



The widths for scalar mesons decay into photon and 
vector meson have the form 



r, 



1 



^1 



M2, 

__ph±\ 

p2 ) 



Mp/^G\ 



X \Ak{p^) M^{ml,p\Ml/J\ 



(44) 



Ml 



27rVp2 



eM„Gx 



X \Bk{p^) ^{ml,p\Ml) 



(45) 



1 



27rvV 



1- 



The expressions (04]), |[45]) and ige]) follow from 



Table 3. Decay data 



,»(1020)-^-,/o 

-r0(lO2O)^^aO 

-r,^(1020)^TJo 
-r,/,(1020)^',a(, 

r 

-1 ao-»77 
-''/o -►77 



(4.40 ±0.21) X 10" 
(7.6 ±0.6) X 10"^ 
6.1 ±0.6 
0.30 ±0.10 keV 
0.3ltn? keV 



m 
m 
m 

[32] 

m 

[32] 



r, 



fo- 



/o^77 



34.2t?|^ MeV 
0.205l!5;2''2 keV 



[3T] 
[31] 



where scalar meson mass is Ms , vector meson mass - My , 
and matrix element ^A is given in l|25p . l(26|) and (|27l) 
respectively. 

Let us discuss difficulties one faces when trying to use 
Eqs. ((36l) - l(4T1) for fixing the couplings Cd,Cm,Cd,Cm and 
mixing angle 9. It is clear that accuracy and even exis- 
tence of relevant experimental data are very important. 
The particle properties are presented in Tableland the 
known decay widths in Table [3] (the latter contains also 
very recent data from KEK [31], however the errors are 
still too big). 

At present there is a big ambiguity in the mass of 
a = /o(600) meson |[32| (although one notices smaller er- 
rors in reference from CLEO [33] in the 4th line in Tabled] 
given for overview purpose only). For that reason we de- 
cided not to use information on cr = /o(600) neither in the 
coupling estimation nor in the width prediction. 

Strictly speaking (|34l) and fSS]) are valid if invariant 
mass of scalar meson is bigger than the pair pro- 

duction threshold 2m k- In fact the physical mass lies 
below threshold, i.e. 2mK > Mf„, Mao (Table [2]). This 
makes impossible to use (|34l) and l(35|) directly in the fit- 
ting procedure (in principle, one can use an approximate 
method [34)). 

From Table [3] one sees that the precision of estimate 
(46) Jqj. Pag^TTTi depends on accuracy with which the total ao 
width Fao^tot is known, as only the ratio 
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is measured [32] ■ Unfortunately Faa^tot has big experimen- 
tal error and therefore extraction of raa-^-nn fron 



leads to a big error. Thus this information should not be 
used in the analysis. 

Further, formally l(36|) and l(40|) are not independent 
as they are expressed through the same factor Ak- We 
prefer to use l(36|l because of a non-trivial assumption ifTTl 
for the couplings used in derivation of (|40l) . For realistic 
values of Fy and Gy the relation ifTTj) is satisfied only 
approximately (see Appendix [Aj . 

From the above reasonings it becomes clear that fixing 
five parameters in question is not an easy task. Therefore 
we reduce the number of independent parameters from five 
to three by applying the large- A^c relations (jB.2p . These 
relations are briefiy discussed in Appendix [BJ Then to find 
values of c^, c„i, 9 one can use l(36|) . l(37|) and (|4T|) . 

In the analysis below we take masses of ao(980) and 



/o(980) equal and put 



980 MeV. Let Ak, Bk 



and i?7r be our estimates for AxiM"^^), BK{M'j^) and 
B^{MjJ. Applying the constraint (|B.2p to ^ we find 
that the scalar mixing angle 9 satisfies the equation 



cos 6* -I- V2 sin 61 = V6 



Bk^ 

A 



K 



where /i = ±1 stands for two possible choices of sign 
in (|B.2p . For the coupHng constants 



Cd 



V2 



K 



-(mlAK - m\RB^), 



((/ - 2ml)AK 
-{p' -2mj,)RB^), (50) 



„2 r,2 



where 



A 



K 



The values of Ak, Bk and -B^ in (|T2l) . extracted 
from l(36|). {37]), (jUJ and experiment, are 



(52) 



decay to 77 or 7 p/w. Secondly, the masses of uq and /o 
mesons are not accurately estabHshed yet. New experi- 
mental results may alter considerably the existing values 
and it is important to know how results of an approach 
depend on the masses. 

Notice that our definition of loop integral ^ auto- 
matically includes the loop kinematic factors, i.e. ^ — 
{a — b)I{a, b) in terms of analytic approach [22]. The quan- 
tity I'f'p is convenient in analysis of loop contribution to 
decay probabilities (UHl), (HII) l(23l) and (|25l) - l[27|) . 

Fig. [4] shows dependence of jif'p on the scalar- meson 
invariant mass. This figure does not include any possible 
interference effects between the pion and kaon loops. From 
Fig. [4] (c) one sees that the pion contribution to fo 77 
decay turns out as large as the kaon one. The loop inte- 
grals crucially depend on the pseudoscalar threshold and 
the relation between the masses of pseudoscalar and vec- 
tor particles. Especially dependencies for the kaon loops 
are complex due to proximity of the KK threshold. Two 
dotted vertical lines in Fig. H] show the physical masses of 
scalar mesons. It is seen that the kaon-loop contribution 
rapidly changes near the KK threshold in the vicinity of 
(49) Mao i-^fo ) ) therefore an error in the mass value may cause 
drastic changes in the KK contribution. 

Note that earlier similar loop integrals for (f) ^ 7 ao //o 
andflo/Zo ^ 7 p/w decays were analyzed in [5,19,22j.In 
particular, authors of pj] concluded that the pion loops 
gave negligible contribution to the decays. 

In this connection we stress that for any observable not 
only the loop integrals but also the coupHngs matter. Thus 
it is important to compare the pion and kaon contributions 
taking into account the corresponding coupHng constants 
as well as the interference effects. Fig.[5]shows the ratio of 
the KK contribution to the total KK + mr contribution 
calculated from ifTTI) . fM)) and Table [H This ratio exhibits 
the effect of interference between pion and kaon loops and 
depends on the ratio Bk/B^^. Our estimates for the cou- 
plings lead to Bk/Btt — 3.37, and arrows in Fig. [5] mark 
the values of relative kaon contribution 



(51) 



Bk w 3.4716 x 10^ MeV^ 
B^ w (1.029 or 4.96) x 10^ MeV^ 
Ak « 2.2456 x 10^ MeV^ 



(here Im B-^ = is assumed) . Inserting these values in l|49p 
and l(50l) one obtains couplings and mixing angle in Ta- 
ble [H The relation between the mixing angle 9 for fi = +1 
and fi — —1 is discussed in Appendix [Bl 



3.2 Analysis of loop integrals 

Here we analyze dependence of the loop inte- 
grals (HI), {Ml, 128} and ([301 on the scalar meson 
invariant mass. Firstly, such dependencies are important 
for processes, where the off-shell scalar resonances enter, 
i.e. for any scalar meson production and its subsequent 



Bk ^im'K,Mj,M'p) 



2 B^ ^{ml, Mj,M^) + Bk if (to^^, M^, M^) 



Bk ^{ml,MlO) 



B^ ^iml,Ml 0) + Bk ^{mj,, M], 0) 



= 1.635, 



= 1.915 



to decays /o 7P and fo 77 respectively. 

Results in Fig. [4] (a) for the decays strongly favor 
the kaon loops compared to the pion loops. That would be 
an argument, additional to OZI suppression rule, for not 
using pion loops for cj) decays. For other processes Fig. [5] 
gives an adequate measure of the pion-kaon concurrence. 
In /o 7 p decay omitting the pion loops would lead to 
« 60% overestimate of the width. The pion loops are very 
important in the two-photon decay of fo: they reduce the 
decay rate by a factor of 1.9. 
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Table 4. Chiral couplings and mixing angle 



B.^, (10'' MeV^) 



1.029 



1.029 



4.96 



.96 



Cd (MeV) 
Cm (MeV) 
(for fi = +1) 
(for = -1) 



-6.39 -52.57 

-58.83 -13.15 

-7.329° 46.271° 

-172.671° 133.729° 



-41.67 -264.37 

-23.93 196.37 

-7.329° 46.271° 

-172.671° 133.729° 




920 940 960 9B0 lODO 1020 1040 



<t> ^ (KK ^) 7 ao 
(a) c^^(KK_^)^fo 
(j) (tt+tt 7 fo 




940 


960 980 


1000 1020 1040 


920 940 960 


980 100 


1020 














ao 


-> {KK 


7 p 








ao 


{KK 


7 


ao — 


* {KK - 


77 


(b) fo 


{KK 


^) 7 p 


(c) fo - 


> {KK - 


77 


fo 


{n'^n^ 


^) 7 p 


fo ^ 


{tt'^tt^ - 


77 


fo 


{KK 


7 LJ 









Figure 4. Comparison of loop integrals squared, ■ Solid line is drawn for kaon loop, dashed line for pion loop. Vertical 
dotted lines mark assumed physical values of scalar meson mass {Mfg and Mag). Two solid lines in (b) account for different 
masses of p and u) mesons. Interference between kaon and pion loops is not included. 




Figure 5. Relative kaon loop contribution, 

\{KK)/{tv+-k' + KK)\'^, to nil) and ^ versus ratio 
Bk/Btt. Curves (b) and (c) correspond to decays fo 7P 
and fo — > 77 respectively (see also legend in Fig. [4|). Our fit 
gives Bk/Btt = 3.37. 



3.3 Model predictions 

In the present model two types of predictions are obtained. 
On the one hand, the observables depend on values of the 



model parameters, and thus can be evaluated after specific 
values are chosen. 

On the other hand, several ratios of the widths turn out 
to be independent of specific values of couplings Cm, Cd, 
Cm, Cd and 9. We find three such ratios which are constant 
in the present model for any values of these parameters: 



r 

-"-an—* 



Co— ♦77 



0.422 



(53) 



GlMaMl-MD 



From experiment (Table[3]) one obtains about 0.93 for this 
ratio. 

Another ratio 



r 

-t aa- 



np 



= 1.043 



(54) 



MM 



{Ml - Ml)Ml 



'P{m 



^ml-Ml-Ml) 



has not been measured so far, though theoretical predic- 
tions exist. In particular, it was shown [l9j that the quark- 
loop mechanism in the two-quark model gives value about 
1/9, the four-quark structure leads to « 0, while kaon-loop 
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mechanism produces almost equal oq ^ 7 p and ao ^ 7 
widths. Our result appears close to the latter prediction. 
From (35), (39) and (43) it is also possible to derive 



r. 



12 



Ml 



(55) 



^2 ,M2 



M2 

ao 



M2 




tf'(m|,M2^,M2) 



Analogous ratio can be deduced from results of [16]. It 
appears to be approximately 5.6. The difference may be 
addressed to dissimilarity between the models, particular 
choice of mass and coupling constant values. 

Now, let us focus on the coupling-dependent results 
calculated according to ([Ml) -{Ml), ^ and (gl]). Tables [5] 
and [6] show predictions of the model and comparison with 
available data. Table [5] presents obtained values of Cd, c,n 
and 9, and strong decay widths for ao and /o mesons. 

In Table [6] one can see results for radiative decay 
widths. The column I shows calculations with our param- 
eters Cd = -6.38 MeV, c™ = -58.83 MeV and = -7.33°. 
The entries in Table [6l which were taken as input in the 
fitting procedure, are marked with asterisk {F^^^fg/F^, 
Fag-,.yy and Ff„^j^). The column II is calculated with 
an "ideal" mixing angle 9 = —35.3° (sm9 = —1/^/3, 
cos 9 = \/2/3) as chosen in [S^. In this case decay a — 
/o(600) TTTT is forbidderlfl, though in fact it should be 
super-allowed. The column III deals with another "ideal" 
mixing angle 9 = 54.7° (0086* = l/VS, sin 6* = ^/V^). 
In this case decay /o(980) — > tttt turns out forbidden in 
contradiction with experimental evidence. Therefore the 
choices II and III do not look reaHstic. 

For comparison of our results with predictions of other 
models we add columns IV, Va and Vb. In partic- 
ular, kaon-loop model [H] (KLM) is selected (column 
IV), which is somewhat similar to the present calcula- 
tion. In columns Va and Vb predictions of vector-meson- 
dominance (VMD) model [21] (Table I therein) are shown 
for two different sets of parameters. The authors apply 
a chiral Lagrangian with strong trilinear scalar-vector- 
vector interaction. 



3.4 Discussion 

In this subsection we briefly compare our results with 
those of KLM [T9] and VMD [24] models and comment 
on correspondence of predicted widths to experiment. 

As it is seen from TableEl our model, contrary to KLM, 
gives not only the rate for decay of a given type {S 77, 
S and (p groups) but also different decay 

rates for ao and /o mesons. For some of the channels the 
results of KLM are qualitative estimates, corrections [14] 



to which should be calculated as discussed in [19]. Nev- 
ertheless, our results for S" ^ 77 and <f) ^ are in 
agreement with KLM within an order of magnitude. For 
the ratio Fag-,.yy/F^-^jao we also get a close value, which 
approximately corresponds to experimental result. At the 
same time we obtain the widths for S jV decays which 
are bigger than the values in the KLM. The latter discrep- 
ancy is due to SU{3) relations for strong interaction (see 
(|A.lip ) and our couplings of (j> and p/uj mesons to KK 
turn out to be different from those used in [19] . 

Regarding the VMD model [24|, one can see from Ta- 
ble [6] that quite a big decay widths for S —f 'yiv are ob- 
tained there compared to our results, while S ^ jp pre- 
dictions differ not so much. Note also a big difference in 
the values of the ratio F^-^jfg / F^^jao ■ 

From results presented in Tables [5] and [6] one concludes 
that predictions for scalar meson decay widths are very 
sensitive to model details. Therefore the future experi- 
ments in which these processes will be studied may help 
to discriminate between different models of scalar mesons. 

In general, for comparison with experiment, in which 
scalar resonances contribute, a more appropriate observ- 
able is the invariant mass distribution. As an example con- 
sider the reaction e+e~ ^ 7* ^ rnr^ at the CM energy 
close to the (/)(1020) mass [2]. This reaction allows for ex- 
traction of the branching ratio 



dB, 



1 d-Z^^—j-TTTT-'^ 



<f). tot 



(56) 



where is the two-pion invariant mass squared. Within 
the present framework this branching ratio can be calcu- 
lated from 



dF 



1. 



dp2 



F^^^f„{p')Bf,^^,{p') X ( - -)lmi?/„(p2). 

(57) 

Here F^^^fg{p'^) is the 4> ~^ ifo decay width (jUJ for 
arbitrary p"^, Dj^ij?) = [p^ -171^ + imfFf„^tot{p'^)]~^ 
is the scalar-meson propagator and branching ratio 
Bfo^TVTrip'^) = Ffg^^^ip"^) / Ffg^totip"^) relates /o -> tttt 
decay width (|32l) to the total fo width Ffg^totip^)- A 
more advanced form of the propagator including both real 
and imaginary parts of the self-energy was suggested re- 
cently [36j. The problem of finite resonance width effects 
in invariant mass distributions for 7r°7r° and tt'^t] in the (j> 
radiative decays is important [37] . 

Note also that in Ref. [19] (Appendix) a more general 
distribution over invariant masses of both the initial and 
final resonances for S jV decays is discussed. 

The scalar octet and singlet mixing angle 9 appeared 
to be crucial parameter in the fit. We should remark that 
a detailed study of mixing angle was performed in [38] 
using the inverse amplitude method. The basic processes 
there were elastic tttt, ttt], KK and Kij scattering and an- 
gle value was different from our estimat^. There are also 
models in which fo is mainly the singlet state with « 0, 
for example, an application of the Bethe-Salpeter equa- 



^ In qq quark model this case corresponds to /o(600) — ss 



* The definition of mixing angle in [38] is also different. 
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Table 5. Strong decays of scalar mesons 





Our estimate (I) 


Estimate (II) 


Estimate (III) 


Exp. value [32] 


Cd, 


MeV 


-6.38 


±32 


±32 




Cm, 5 


MeV 


-58.83 


±42 


±42 






6 


-7.33° 


-35.26° 


54.73° 




r 

-t aQ—>7rri-, 

-^/o— 

raQ.tot 5 


MeV 
MeV 
MeV 


14.2 
41.8 
17.8 


172.4 
775.7 
215.4 


172.4 

7.4 X 10"^ 
215.4 


•AA 0+22.7 
50 - 100 



Table 6. Model predictions and available data for electromagnetic decays. (Asterisk * marks the experimental values, used for 
extraction of couplings. For couplings used in columns (I), (II) and (III) see Tabled 



Observable 


Our est. (I) 


Est. (II) 


Est. (Ill) 


m (IV) 


m (Va) 


[24] (Vb) 


Exp. value [32] 


EpnjiSL 10-4 

^ <f>,tot 


1.67 


2.13 


2.13 


1.4 






(7.6 ±0.6) X 10"^ 


£p^2lSL 10-4 


4.40* 


2.31 


4.63 


1.4 


4.92 ±0.07 


4.92 ±0.07 


4.40 ±0.21 


-T^— 700 


2.64 


1.08 


2.17 


1 


0.26 ±0.06 


0.46 ±0.09 


6.1 ±0.6 


Fo,Q — ^77 5 keV 
-Oo^77i keV 


0.30* 
0.31* 


0.383 
0.323 


0.383 
0.62 


0.24 
0.24 


0.28 ± 0.09* 
0.39 ±0.13* 


0.28 ±0.09* 
0.39 ±0.13* 


0.30 ± 10 
31+0 08 

'-'•'^^-0.07 


Fo,Q — ^7P7 keV 
rfo^jp, keV 

-^ao — >^LO J 


9.1 
9.6 
8.7 
15.0 


11.65 
0.95 
11.15 
7.93 


11.65 
16.6 
11.15 
15.85 


3.4 
3.4 
3.4 
3.4 


3.0 ± 1.0 

19 ±5 
641 ± 87 
126 ± 20 


3.0 ± 1.0 
3.3 ±2.0 
641 ± 87 
88± 17 





tion with a linear confinement qq potential to calculation 
of scalar- meson mass spectrum [39]. 

This discrepancies, to our opinion, may not be caused 
only by differences in the appHed models. The problem 
may be related to a non-trivial structure and behavior of 
light scalar resonances. As it was emphasized by Bugg 0, 
unification of observed resonances in elastic scattering ex- 
periments with corresponding ones seen in radiative de- 
cays is a big challenge. They show up in different ways, 
and it would be important to build a consistent bridge 
between these and those properties of resonances. 

3.5 Possible interactions beyond the model 

Comparison of our results with predictions of other mod- 
els, especially the results independent of the choice of 
couplings, shows that the present model does not allow 
one to reproduce the ratio Fao^jp/ Fao^jLj — 1/9, which 
is obtained in qq model and quark-loop mechanism [19j, 
or Faa^-yp/ Fag^-yu ~ in qqqq model ^9j. The present 
model is insensitive to the structure of scalars. 

At this point one can think of a direct (or contact) 
coupling of scalars to two photons as an extension of the 
present model. Similar terms were introduced in [35] and 
have the order O(p^) 

Li - g (5°^*/r/+M-) + (/r/+M-) - (58) 

where g, g' are coupling constants and f!f'^ is defined in 



Analogously, constructing the C and P invariant terms 
with y^i/ one can propose the SjV interactions 

L2 = g" (5°^*/r V) + 5"'^""^ (f+V^u) (59) 

which are bilinear in resonance fields. Lagrangian l(59|) has 
the order C'(p^) and contains two more couplings g", g'" . 
Of course, these terms do not violate chiral symmetry. 
Four additional coupling constants g,g',g",g"' should be 
fixed from certain observables. 

To our opinion, Lagrangians (|58|) and l|59p can be use- 
ful in phenomenological descriptions of scalar radiative de- 
cays. They may represent effects related to specific quark 
structure of the scalar mesons, which is not accounted for 
in chiral Lagrangian IjB.ip . This aspect Hes beyond the 
scope of the present paper. 

4 Conclusions 

Within ChPT with vector and scalar mesons [25] we 
have calculated the radiative decays ao ^ 77, /o ^ 77, 
(j> 700 and (j> 7/0. These decays and corresponding 
invariant mass distributions can be measured in e+e~ an- 
nihilation in Frascati by KLOE [2T1I40] and Novosibirsk 
with VEPP-2000. 

The derivative and non-derivative couplings of scalar 
mesons to pseudoscalar ones are consistently included. 
The gauge invariance of the amplitudes and cancela- 
tion of divergencies from different diagrams are explicitly 
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demonstrated. The obtained ampHtudes are finite with- 
out counter terms. For (/> decays, in addition, we used the 
relation Fy — 2 Gy between electromagnetic and strong 
couplings of vector mesons in order to get rid of divergen- 
cies. This relation was previously discussed in ^30j in con- 
nection with alternative approaches. Hidden Local Gauge 
Symmetry [28] and massive Yang-Mills models [29]. Note, 
that this relation does not follow from chiral symmetry 
but does not contradict it as well |30j . 

The scalar flavor singlet-octet mixing angle is ob- 
tained from the fit, as well as estimates for octet chiral 
couplings Cm and Cd- It should be noted that the values of 
these parameters strongly correlate with the mixing angle. 

For the flavor singlet couplings Cd.m we relied on the 
relations Cd^m = Cd,m/V^ to flavor octet couplings in the 
large- A'^c limit [25]. One may argue whether the large- A^c 
consideration is applicable to scalar mesons, especially in 
view of Unitarized ChPT results [18]. In this connection, 
a flt without any large- A'c restriction would be an exten- 
sion of the present approach. However difficulties related 
to the fitting procedure may arise, in particular, two more 
free parameters Cd,m appear, and in view of scarce exper- 
imental data a non-trivial procedure is needed to reduce 
ambiguities in the results. 

In the present model we obtained the widths of ao(980) 
and /o(980) decays: ra„,tot = 17.8 MeV, /"qo^ttt? = 
14.2 MeV, r/o_^^^ = 41.8 MeV. Many of the calculated 
observables are in satisfactory agreement with experiment. 
At the same time the calculated ratio /^-.^/q //0_,^ao — 
2.64 only qualitatively agrees with the experimental value 
6.1. The results of the present approach are also compared 
with those of previously developed kaon-loop model [19] 
and vector- meson-dominance model [24]. 

Predictions for the widths of ao(980) and /o(980) de- 
cays into 7p(770) and 7Li;(782) are also given (see Table[6]). 
The processes, to our opinion, are of interest for experi- 
mental programs in Jiilich with COSY [20j and R-ascati 
with DA^NE (or its upgrade) [4T1I42] . 

Within the present model and one-loop approxi- 
mation we found several ratios of the widths which 
are independent of the coupHngs constants. Namely, 
Faa^-^i-y / r^^^ao = 0.422, which is in a qualitative corre- 
spondence with experiment, Fag^^p/ Fag^^uj = 1.043 and 
rao~*^p(uj)/ F^-,^ao ~ 12 which have not been tested ex- 
perimentally yet. 

Our calculations show that many predictions are in 
agreement with experiment, and therefore support the as- 
sumption that ao(980) and /o(980) fit in the Hghtest scalar 
meson nonet. However, it is difficult to make an unambigu- 
ous conclusion. 

The present work makes a solid ground for further 
studies of scalar mesons, not only the lightest ones ao(980) 
and /o(980). The model can be applied in processes of 
two-photon production of hadronic states with intermedi- 
ate scalar resonances. These processes occur in nucleon- 
nucleon and electron-positron collisions (like e+e^ 
e~'"e~7r"'"7r~). 
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A Chiral Lagrangian for pseudoscalar and 
vector mesons 



In calculations we use 0{p^) ChPT Lagrangian for 
pseudoscalar mesons vector mesons and photons, de- 
rived by Ecker et. al. [25], where spin-1 mesons are de- 
scribed by antisymmetric tensor fields V^^. This La- 
grangian has 0(j3^) chiral power in sense of its equiva- 
lence to the ChPT Lagrangian in which no explicit res- 
onances are introduced (see [25,30] for details). In the 
present problem it is sufficient to keep 

+ ^ {V,.fn + ^ ( V^^'O , (A.l) 

where (• • • ) stands for the trace in fiavor space. The pion 
weak decay constant ~ 92.4 MeV, Fy and Gy are cou- 
pling constants. The electromagnetic field B^^ is included 
as an external source, F^i, — dp,Bi, — d^Bp, is the electro- 
magnetic field tensor. The quark mass matrix 

X = 2Bo diag(m„, nid, nis) (A. 2) 

is expressed in terms of light quark (g„, qd,qs) masses and 
chiral condensate: {0\ququ\0) = — f^Bo{l+0{mq)) . In the 
limit of exact isospin symmetry x — diag(mj, mj, 2m\ — 
m^). The terms in Lagrangian relevant for scalar meson 
sector are discussed in Appendix [BJ and are denoted as 
in present paper. 

Pseudoscalar meson nonet {J^ = 0^) contains the 
^flavor octet of Goldstoue bosons and the 1 flavor sin- 
glet, namely 770 field. We combine singlet and octet into 
nonet following f44]. Thus fiavor SU{3) multiplet for pseu- 
doscalar mesons is 

<P = -^(ttiAi -I- ^2X2 + T^sh 
v2 

+K1X4 + K2X5 + K3Xe + K4X7 

+m^s + VoXo), (A. 3) 
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where \a {a — 1, •■•,8) are the Gell-Mann matrices, Aq = 



1, and the physical fields are defined as 



TT^ - -^(^1 T *7r2), = -^(ifi T ^K2), (A.4) 



Of course such a scheme is well-defined only approximately 
if the U{1) axial anomaly is neglected. We do not omit 
problematic ij meson within present approach, as it is in- 
volved in the dominant decay of uq meson ag Trry. For 
T] - rj' mixing we choose the two-parameter scheme [45j : 



rj = cos Vs - sm Oq tjq, 
77' = sm9sT]s + cosdor]o. 



(A.5) 



Note that i] and 77' are not the orthogonal states. The 
angles 6*0 = —9.2° and 6*8 = —21.2° are discussed and 
determined in |45P46j from experiment. 

Further, f^'^ in (|A.ip for the external electromagnetic 
field reads 



where 



exp 



/ 1$ 



(A.6) 



(A.7) 



carries non-linear parametrization of the pseudoscalar 
field. The quark charge matrix is Q = diag(|, — |, — -i) = 
^Aa-l-^^As, and the electron charge is e = ^/ina ~ 0.303. 

The definition oiD^JJ, XI^V^" and u" in ([Al]) can be 
found in the original work |25) . 

From \K.\\ the following interactions for physical 
fields in Oij?) order can be produced [43]: 



£^PP = ~ieB,,{TT+ d,, ^' +K+ df^K^), (A.8) 
Cy^4.f = e^B'^Bf^iir+n- + K+R-), (A.9) 



where for any a and b notation a 9^ 6 = a 9^ 6 — 6 9^ a is 
introduced. 



(A.IO) 



C 



G 



VPP 



z-f [pO,(2a^^+9'^^- + d^'K+d^K-) 

J TT 

+(f>^,{-V2dt'K+d''K-)]. (A.ll) 



<?!)(1020) 



(?i(1020) 



P- 



P+ 



-Gy _ ;+,-] 



(^6(1020) 



p- 
p+ 



p- 



(M) 

X 

('^) p- 



Figure 6. The C'(p^) vertices from chiral Lagrangian l|A.ip 
Dash line stands for kaon, double solid - for vector meson (j) 
wavy line - for photon. 



-yVPP 



eFy 

J 7T 

eFy 
2eGv 



(A.12) 

K+K- 



eGy 

From these terms one can derive the 0{p^) vertex func- 
tions shown in Fig. [6l 

Chiral couplings Fy and Gy can be extracted from 
the vector meson partial widths (see [.25.4.3j ). From l|A.10p 
and IjA.lip one calculates the following decay widths: 



r 

p — *7r7i 



F 



487r/4^ ' 



4m 



2\3/2 



p— ^e^" e 



127rM„ 



(A.13) 
(A.14) 



For the decays (j) KK, uj e+e^ and (f) e+e^ one 
has to take into account the SU (3) relations for strong and 
electromagnetic coupHngs implemented in l|A.8P - (|A.12p . 
Tables [7] and [8] show values for Gy and Fy which are 
obtained from experimental widths. 

The condition Fy = 2 Gy is important for amplitudes 
to converge (see Section 12. 3p . From Tables [7| and [8] one 
sees that this relation is satisfied only approximately, and 



S. A. Ivashyn and A. Yu. Korchin: Radiative decays with light scalar mesons and singlet-octet mixing in ChPT 
Table 7. Values of electromagnetic coupling constants for vector mesons 



15 



rp^^+^-, keV Fv (MeV) 
7.02 ±0.11 156.162 


r^^e+e-, keV Fv (MeV) 
0.60 ± 0.02 137.629 


F^_,+,-,keV Fv (MeV) 
1.27 ±0.04 161.629 


^^P^M+z.-' keV Fv (MeV) 
6.66 ±0.20 152.358 


^..^P+M-'keV (MeV) 
0.76 ±0.26 154.98 


F^_,^+,,-, keV Fv (MeV) 
1.21 ±0.08 157.738 



Table 8. Values of vector-meson coupling to two pseudoscalar 
mesons (all values are in MeV) 



with scalar mesons is defined by 







TT + TT 


K+K- 






exp. width: 


146.4 








Gv- 


65.183 








2 Gv: 


130.366 








exp. width: 


0.144 










(suppressed) 








exp. width: 




2.096 


1.448 




Gv'- 




53.09 


54.45 




2 Gv: 




106.18 


108.9 



the closest values are obtained from (|A.13P and l|A.14p for 

the decays p tt+tt" and p e+e". 



B Chiral Lagrangian for light scalar mesons 



The O(p^) ChPT Lagrangian, which explicitly incor- 
porates scalar mesons and their interactions with pseu- 
doscalars reads [25| 



tB _ jChPT I jChPT 
^ ^ ^oktet -^singlet! 



(B.3) 



rChPT _ '^^d 
^oktet £2 
Jtt 



2ecd 

J TT 

2e' Crf 
P 



-ChPT 

■^singlet P2 
Jtt 



f2 

J TT 

P 



J TT J TT 



(B.l) 



X+ = u'^X'^'^ + uxu. For other notation and definition see 
Appendix [A] and |25) . 

Scalar octet 5°'^* and singlet have a priori inde- 

pendent couplings Cd,c„i and those with hats Cd,Cm- Nu- 
merical values of these couplings are determined by the 
underlying QCD. However it is difficult to find Cd, Cm, 
Cd and Cm at energies about 1 GeV because of the non- 
perturbative regime of QCD. From assumption of large 
number of quark colors {Nc - 
octet and siglet (with "tilde" 
tions 



oo) it was shown [25] that 
chiral couplings obey rela- 



^/3' 



Cd 



Cd 



(B.2) 



Applicability of (|B.2p to scalar meson radiative decays 
gives rise to some doubts ([E], for instance). Anyway we 
use these constraints to reduce the number of independent 
parameters in Section ISTTl 

For description of scalar meson radiative decays we 
expand in (|B.ip in series in ^. The 0{p^) interaction 



Apparently Lagrangian (|B.3P does not yield direct contact 
coupling of scalar meson to two photons. 



In order to apply l|B.3P to the physical scalar fields one 
has to assume certain multiplet decomposition ([1]) . Which 
has to be consistent with phenomenology. The prominent 
feature of ao is its dominant decay to Trrj. The KLOE [1] 
showed almost no contribution of /o(980) resonance in 
comparison with ao(980) in the reaction (j)[\Q2Q) jn^rj. 
Thus we suppose that isovector ao(980) and isoscalar 
/o(980) do not mix with each other. Violation of isospin 
conservation, related to a possible ao-/o mixing, is a sub- 
ject for a separate work. This issue can be studied for 
example by means of the dd {'^He Gq ^) '^He ir'^rj 
reaction at COSY [47]. 

As far as we are interested in physical scalar fields, 
which are combinations of singlet and octet states, it 
is convenient to introduce effective couplings gs -- con- 
structed from constants Cd,Cm,Cd,Cm- This allows one to 
rewrite Lagrangian in a simpler form. 

Let S stand for any scalar field, ao,/o or cr, and P - 
for pseudoscalar n— 7r°,7r='= or K"^, K^, K^. Then La- 
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' 9SPP 
^. Jtt 



_19SPP 



Jtt 



(/Or 



P+ 



(2 
Jt: 



■' ''■e9SPPi,+ ,-^ 

N JTT 



Figure 7. The O(p^) vertices corresponding to La- 
grangian l|B.4p . Dotted line stands for scalar meson S, dashed 
line - for pseudoscalar P. Couplings are shown in Tabled see 
also HB.SII. 



(B.4) 



grangian (|B.3P can be reduced to 



9S7, 



9Srjri 2 I 
+ dSTT-n-^ V 



+gsKK{K+K- + {-lY^K''K°) + {gs^^/2){df, n)^ 

+gsKK{d^,K+d^K- + i-iy^d^.K^d'^K") 
+5S77r7reS^7r+ 5p tt" + 9SfKKeBf^K+ 5^ K~ 
+9S'yyrr^e^B^B''7r+TT- + gs-y■yKKe^B^B^'K+K-^ , 

where Is = for /o, a and Is — ^ for oq. We introduced 
the effective coupHngs gsTi-K, gsriri, etc. Hsted in Table [9l 
The coupHngs which are absent in Table [9] are equal to 
zero. In addition, for any scalar meson S the following 
relations for electromagnetic couplings hold 



gSyirTz — ^gS-TTTz^ 

gSiKK = —igsKK, 

gSf-yTTTT gSTTTT: 

gSj-yKK = gSKK- 



(B.5) 



Note also that 



z = 



cos Oq 



-V2 



sm Og, 



COS(«8 - fo) 

1.53, 



(B.6) 



where denominator in Z is equal to determinant of the 
transition matrix (|A.5P from (ryg,77o) to {rj^rj'). 

Lagrangian l|B.4p leads to the vertices shown in Fig.[7l 
Table [9] shows expressions for the effective couplings 
as well as the corresponding chiral powers. Some of the 
couplings include the masses of Goldstone bosons and are 
C(p^), while the others are 0{p^). Of course each term in 
Lagrangian (|B.4p carries power 0(p^). 

Let us now make a remark on relation between the sign 
of the parameter fj, in (|B.2p and scalar- meson mixing angle 
in ([1]) . As long as present consideration does not involve 
tr meson we can drop the relation for cr in (HJ and observe 



Table 9. Effective couplings and their chiral powers for scalar 
mesons 



ff/TTTT = -m^(4cm cos 61 - 2^2/ \/3 c„i sin 61), 
ff/,,,, = -4/3 Cm(4m|- - m^) 008 61 

-2^/(3^3) c™(8m^ -5m^) sine), 
QfKK = —m|f (4 Cm cos 6* + V2/\/3 Cm sin 6*) . 



P/tttt = 4cdCose — 2\/2/\/3ctjsin6', 
Qfnrj = 4: Cd cos 6 + 2\/2 / Cd sin 9 , 
QfKK = 4:CdCOs9 + \^/\^Cdsm9. 



0{P°) 



g^TVTv = -mKiCm sin 9 + 2y/2/y^ Cm cos 9), 
Qann = -4/3cm(4m|f - m^) sin6' 

+2^/ (3^3) Cm (8mif - Sm^ ) cos 6») , 
gcrKK = -m^ (4 Cm sin 6* - V2/V3 Cm cos 6*). 



gaiT-K = AcdSin9 + 2-^/2 /^/icd COS 9, 

garir, = 4: Cd sin 9 — 2\^ / Cd COS 9 , 

gcrKK = AcdSin9 - ^/2/^/?>CdCos9. 



0(p") 



gaKK 
ga-Kri 



-\/2cmmjf 
-2Z^f2l^ 



QaKK 



y/2cd , 

2ZV2/V3cd 



0{p") 



a nontrivial property: the change ^ — + — ^ is equivalenlH 
to the change 9 ^ ir — 9. 



C Modifications of the model for virtual 
photons 



The complete set of ©(p^) diagrams has to incorporate 
all contributions determined by 0{p'^) Lagrangians IjA.SP - 
(|A.12p and l|B.4p . Lagrangian I|A.1QP generates electro- 
magnetic form factors (FF's) for pseudoscalar particles 
inside the loops. These FF's should replace the tree-level 
PP7 vertices marked by arrows in Fig. [J, if the FF's do 
not increase the chiral power of a diagram and they are 
calculated from 0{p'^) Lagrangian. 

Note that electromagnetic FF's of kaons and pions 
have been studied in various approaches (let us just men- 
tion the considerations for on-mass-shell pions [48lj49j and 
kaons [43j). The FF calculated from ChPT Lagrangian in- 
cludes direct photon-vector meson transition, i.e. vector 
meson dominance, as well as ordinary contact interaction 
(see illustration in Fig. [8] and Appendix [Aj . Fortunately 
the real photons do not couple to vector mesons within this 
approach (see, e.g., [43j for discussion of this and one-loop 
modification of electromagnetic vertex) . Therefore as long 
as one is interested in processes with real photons there 
are no 0{p'^) diagrams additional to those shown in Fig.dJ 
and therefore Fig. [T] gives the complete set of diagrams in 
this order. The similar reasoning applies to consideration 



This equivalence is reflected in Table S) 
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Figure 8. The ©(p^) electromagnetic vertex of (off-mass-shell) pseudoscalar meson in ChPT. All possible intermediate vector 
resonances V = p°, w, 0, ... in general contribute. For real photons only the first term on the r.h.s. is non-zero. 



Table 10. Table of typical D-dimensional integrals (for any 
vector and complex number -R). 



X [J,-ln(2i^)] 

X [/, + l-ln(fii^)] , 

J {27r)-" (i^-2i<3+fl):i 2(4ir)^ Q^-H 

X {0,0,5"", (s-^Q" 



of diagrams shown in Fig. [2] and Fig. [3] for (/) ^ 7ao//o 
and /o/oo ^ 7 p/^ decay respectively. 



D Dimensional regularization. Loop integrals 



In calculation of loop integrals we use the dimensional 
regularization method (see, for instance, § 7 and Ap- 
pendix [B] in pj). 

The dimension of space-time D — 4 — 2e in the limit 
e corresponds to that of 4-dimensional Minkowsky 
space. In the text this limit is assumed in all expressions. 
Integration measure for 4-dimensional space is replaced 
by that for D-dimensional space: d*q {A^Yd^q, where 
arbitrary regularization parameter A has units of mass. 
Integrals with this measure are defined via the analytical 
continuation from the space with the integer number of di- 
mensions. Metric tensor obeys the condition g^'^g^u = D. 
The Dirac matrices satisfy the anti-commutation relations 
{l^il'^} — "^g^^ ■ Here 7'^7p — D, and the ordinary trace 
formulae are generalized to Tr(7^7'') = 2^/'^gt"', etc. 

In calculation of loop diagrams the typical integrals 
presented in Table [TOl arise. For divergent terms we define 

/e = 1/e - 7e + ln47r. 

Euler-Mascheroni constant 7e ~ 0.57721566490 can be 
expressed in terms of gamma-function derivative 7c = 

—r'{l)—~J^ (ia;exp (— x) Ina;. 
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